Quantum gases in trimerized kagome lattices 
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We study low temperature properties of atomic gases in trimerized optical kagome lattices. The 
laser arrangements that can be used to create these lattices are briefly described. We also present 
explicit results for the coupling constants of the generalized Hubbard models that can be realized in 
such lattices. In the case of a single component Bose gas the existence of a Mott insulator phase with 
fractional numbers of particles per trimer is verified in a mean field approach. The main emphasis 
of the paper is on an atomic spinless interacting Fermi gas in the trimerized kagome lattice with 
two fermions per site. This system is shown to be described by a quantum spin 1/2 model on the 
triangular lattice with couplings that depend on the bond directions. We investigate this model by 
means of exact diagonalization. Our key finding is that the system exhibits non-standard properties 
of a quantum spin-liquid crystal: it combines planar antiferromagnetic order in the ground state 
with an exceptionally large number of low energy excitations. The possibilities of experimental 
verification of our theoretical results are critically discussed. 



PACS numbers: 03.75.Fi,05.30. Jp 

I. INTRODUCTION 

The experimental realization of the Bose-Einstein con- 
densation (BEC) [1] linked the physics of cold atoms with 
that of weakly interacting many-body systems, tradition- 
ally studied by condensed matter physics. More recently, 
the seminal theory paper by Jaksch et al. [2] , followed by 
equally seminal experiments by Greiner et al. [3] on the 
Mott insulator (MI) to superfluid (SF) transition have 
paved the way towards the analysis of strongly corre- 
lated systems within the physics of cold atoms. In this 
sense, the physics of cold atoms is nowadays merging 
with condensed matter physics, solid state physics, and 
quantum information at the same common frontiers and 
open challenging problems, such as for instance BEC- 
BCS crossover (see for instance [4]), fractional quantum 
Hall effect (c.f. [5]), physics of ID systems [6], etc. Quan- 
tum information has given new impulses towards the un- 
derstanding of quantum phase transitions [7] , and to un- 
derstand better the known (and develop new) numerical 
methods of treating many-body systems [8] . 



A. Atomic lattice gases 

One of the most fascinating playgrounds of the cold 
atom physics is provided by ultra cold lattice gases, i.e. 
cold atoms trapped in optical lattices produced by stand- 
ing laser waves where, in the case of red (blue) detuned 
lasers light, the potential minima coincide with the in- 
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tensity maxima (minima) [9] . This technique has been of 
enormous interest during the past years. The setup can 
be chosen to be one-, two- or three-dimensional, where 
the lattices form ranges from simple periodic (such as a 
square in 2D, or cubic lattice in 3D, respectively) to more 
exotic lattices, such as hexagonal [10], or kagome [11] lat- 
tices, created with the use of superlattice techniques [12]. 
In experiments, optical lattices offer an unprecedentedly 
wide range of tunable parameters, which can be changed 
during the evolution in situ and "m vivo" , i.e. in real 
time. These possibilities, on one hand, link strongly the 
physics of ultracold atoms in optical lattices to various 
areas of condensed matter physics, and on the other, they 
open completely new ways to study quantum many-body 
systems, to perform in various ways quantum information 
processing (cf. [13, 14]), and even to realize special pur- 
pose quantum computers, so-called quantum simulators 
[15]. 

The physics of ultracold atomic gases in optical lattices 
is in general described by various versions of the Hub- 
bard model, which is probably the most important and 
structurally simple model of condensed matter physics, 
capable nevertheless to describe an enormous variety of 
physical phenomena and effects [16, 17]. Atomic ultra- 
cold gases may serve as a "Hubbard model tool- kit" [18] , 
and several models have been discussed in more detail in 
this context: the most simple Bose Hubbard model [2] 
(for the seminal condensed matter treatment see [19]), 
the Fermi-Fermi model (which should eventually allow for 
quantum simulations of high T c superconductivity [20]) 
the Fermi-Bose model (which leads to creation of compos- 
ite fermions via fermion-boson, or fermion-bosonic hole 
pairing, c.f. [21, 22]), or Bose-Bose, or more generally 
multicomponent systems. Quenched disorder may be in- 
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troduccd in a controlled way to such systems [23, 24], 
which opens the possibility of studying the physics of 
disordered systems in this framework. 

In a certain limit Hubbard models reduce to spin mod- 
els and this possibility has been also intensively investi- 
gated recently for both atomic gases [10, 25, 26], and 
ion chains [27]. Spin models enjoy particular interest be- 
cause of their simplicity and thus possible applicability 
in quantum information processing (c.f. [26, 28]). In this 
paper we will discuss yet another possibility, i.e. the pos- 
sibility of studying frustrated quantum antiferromagnets 
(AFM). 



B. Quantum antiferromagnets 

Quantum antiferromagnets, and in particular frus- 
trated AFM's are in the center of interest of modern con- 
densed matter physics (for a recent review, see [29]). One 
of the reasons for this interest is that frustrated AFM's 
are believed to explain certain aspects of the high T c su- 
perconductivity [16]. In this context frustrated spin- 1/2 
models have attracted particular attention. At the same 
time almost all of these models are notoriously difficult to 
handle analytically and numerically. The only exceptions 
are those models that exhibit long range Neel type order, 
since there is a powerful method by which long range 
order can be identified numerically (see e.g. [30]), and 
if it exists, the semiclassical (spin-wave) approximation 
yields satisfactory results. In 2D only very few exactly 
solvable spin-1/2 models are known [29]. In ID exact 
results can be obtained by the Bethe-Ansatz technique 
in a number of cases [31]. Moreover non-perturbative 
bosonisation techniques and powerful numerical meth- 
ods such as the Density Matrix Renormalization Group 
(DMRG) method can be applied. However, DMRG tech- 
niques become very difficult to handle in the case of dis- 
ordered systems [32]. In 2D, in the absence of long-range 
order, apart from renormalization group approaches nu- 
merical methods offer the only possibility to investigate 
frustrated spin systems. However, Quantum Monte Carlo 
(QMC) simulations of Hcisenberg AFMs on frustrated 
lattices, such as the triangular and the kagome lattice, 
suffer from the "negative sign" problem. For instance, 
attempts to obtain useful results for the Hcisenberg AF 
on a triangular lattice (TAF) by QMC have been futile 
as a consequence of this problem. Because of the expe- 
rience with this and other frustrated models, we expect 
that the "negative sign" problem also invalidates QMC 
for the system to be studied in this paper. In contrast 
with the failure of QMC, exact diagonalization of the 
Hamiltonian of the TAF for rather small cells of the lat- 
tice has produced the main result for this model: its 
ground state, contrary to earlier conjectures, was shown 
to be long range ordered (see [33]). 

According to C. Lhuillier and her collaborators quan- 
tum Hcisenberg AFM's at very low temperatures exhibit 
4 distinct kinds of quantum phases: 



• Semiclassical ordered Neel phases, characterized by 
long range order in spin-spin correlation function, 
breaking of the SU{2) symmetry, and gapless spec- 
trum with AS Z = 1 magnon excitations. The 
standard example of such order is provided by the 
Heisenberg AFM on a square lattice in 2D. The the- 
oretical description of such systems using the spin 
wave theory (cf. [16]) is quite accurate. 

• Valence Bond Crystals (VBC) (or Solids), charac- 
terized by long range order in dimer coverings, with 
prominent examples being the AKLT model in ID 
[34], or the Heisenberg model on the 2D checker- 
board lattice [29, 35] (corresponding to a 2D slice 
of a pyrochlore lattice). VBC's exhibit no SU{2) 
symmetry breaking, short range spin-spin correla- 
tions, long range dimer-dimcr order and/or order 
or long range order of larger S = plaqucttes, and 
gapped excitations in all S sectors. 

• Resonating Valence Bond spin liquids (Type I), 
that exhibit a unique ground state, no symme- 
try breaking of any kind, gapped fractionalizcd 
"spinon" excitations, and vanishing correlations in 
any local order parameter. An example of such 
a spin liquid is realized in the, so called, ring ex- 
change model on the triangular lattice [29]. 

• Resonating Valence Bond spin liquids (Type II), 
that exhibit no symmetry breaking, no long range 
correlations in any local order parameter, and an 
extraordinary density of states in each total S sec- 
tor. Numerical work by Dommange et al. [36] sup- 
ports the conjecture of gapless deconfined "spinon" 
excitations in this scenario. An example of such a 
spin liquid is believed to be realized by the Heisen- 
berg spin 1/2 model on the kagome lattice [29, 37- 
43]. 

The kagome spin 1/2 antiferromagnet (KAF) seems to 
be a paradigmatic example of type II RVB spin liquids, 
but unfortunately so far no experimental realization of 
this model has been found among solid state systems. 
Only the spin 1 KAF can be realized in solid state ex- 
periments, but that system has a gap to all excitations, 
i.e. it does not belong the type II spin liquids [44]. The 
physics of the spin 1/2 KAF is, however, not yet fully 
understood. There are papers that suggest VBC type 
order with large unit cells [45]. 

C. Spinless interacting Fermi gas in a kagome 
lattice 

We have proposed recently how two realize the trimer- 
ized kagome optical lattices using superlattice techniques, 
and have studied various kinds of quantum gases in such 
lattices [11]: i) a single component (polarized) Bose gas, 

ii) a single component (polarized) interacting Fermi gas, 

iii) a two component ("spin" 1/2) Fermi- Fermi mixture. 
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In the subsequent paper [46] we have concentrated on the 
second of the above mentioned situations and studied the 
polarized interacting Fermi gas in the trimerized kagomc 
lattice at the filling v = 2/3. Using the method of ex- 
act diagonalization of the Hamiltonian we have shown 
that the system exhibits novel kind of behavior at low 
temperatures, which has led us to propose a new class of 
possible behavior of frustrated AFM's: 

• Quantum spin-liquid crystal, characterized by the 
long range Neel type of ordering at low T, gapless 
spectrum, and anomalously large density of low en- 
ergy excitations. 

This paper is devoted to the presentation of the de- 
tails of the theory described in above mentioned two let- 
ters, Refs. [11, 46]. First, we discuss briefly the general 
properties of interactions in trimerized kagome lattices 
as well as the case of a single component Bose gas in the 
trimerized kagome lattice. Then we focus, however, our 
attention on a trimerized kagome lattice loaded with a 
spinless Fermi gas with nearest-neighbor interaction. At 
2/3 filling per trimer such Fermi gas behaves as a frus- 
trated quantum anti-ferromagnet, and exhibits quantum 
spin-liquid crystal behavior. The motivation to study 
this model is at least fourfold: 

i) In a magnetic field such that the trimerized KAF is 
driven into the magnetization plateau at 1/3 of the satu- 
ration magnetization, the physics of the KAF is described 
precisely by our model [37-39] . Studying our model will 
thus exactly shed light on the theory of KAF and, hope- 
fully, on experiments on the KAF. 

ii) Theoretical studies (using exact diagonalization of 
the Hamiltonian) indicate that the model has the fasci- 
nating properties of, what we have termed a quantum 
spin-liquid crystal. We expect the behavior observed 
in this system indeed to be generic for other "multi"- 
merized systems. First of all it is clear that optical meth- 
ods allow for creating many similar spin models with cou- 
plings depending on bond directions. In the simplest case 
this can be accomplished for a square lattice where one 
could achieve a "square lattice of small squares" , for the 
triangular lattice to obtain a "triangular lattice of small 
triangles" etc. One can expect that when such proce- 
dures are realized for frustrated systems, this might lead 
to similar effects as for the kagome lattice. 

hi) One of the most fascinating possibilities provided 
by the optical lattices is the possibility of "on line" mod- 
ifications of the lattice geometry. We may go from trian- 
gular to kagomc lattice in real time in a controlled way. 
Trimerization (or generally "multi"-merization) is a new 
experimental option, and it is highly desirable to explore 
its consequences. Our model (apart from the model of 
the Bose gas in the trimerized kagome lattice) is one of 
the simplest ones to explore these consequences. 

iv) Last, but not least the model is experimentally fea- 
sible. 



D. Structure of the paper 

The paper is organized as follows. In section II we 
briefly describe the laser arrangement that can be used 
to create a trimerized optical kagome lattice. In subsec- 
tion II B we first introduce the Hamiltonian that governs 
the particle dynamics in the lattice: it is a generalized 
Hubbard model that can be used as a model for bosons, 
fermions, as well as for boson-fermion, or fermion-fermion 
mixtures in the lattice. We show under which conditions 
a tight-binding description of the particle dynamics is ap- 
propriate in such a lattice, and present results of the cal- 
culations of the Hubbard model couplings as a function 
of parameters of the systems and the degree of trimer- 
ization. In the next subsection (II C) we present in some 
detail the results concerning the physics of a Bose gas in 
the trimerized kagomc lattice. Here we generalize the re- 
sults of Ref. [11], obtained in the hard core boson limit, 
to the case when more than one boson can be present 
at the same lattice site. The last subsection (II D) dis- 
cusses in short the case of a Fermi-Fermi mixture in the 
trimerized kagome lattice. 

In section III we start our discussion of the case of 
2/3 filling of the trimerized kagome lattice with spinless 
fermions. We focus our attention on the case of strong 
intra trimer and weak inter trimer coupling. First, we 
discuss various methods of creating an ultracold polar- 
ized interacting Fermi gas in an optical lattice (subsec- 
tion III A). Then we discuss in detail the intra-trimcr 
dynamics. In the following subsection (subsection IIIB) 
we show that the low energy physics of such a gas at 2/3 
filling is described by an effective spin-1/2 Hamiltonian 
with strongly anisotropic couplings. The exchange con- 
stant J of this Hamiltonian is proportional to the inter 
trimer atomic interaction potential which, in the low- 
energy limit, can be attractive or repulsive, depending 
on the species of interacting atoms. In favorable cases it 
can also be manipulated by a magnetic Fcshbach reso- 
nance. The relation of the model to the Heisenberg spin 
1/2 AFM in the kagome lattice is discussed in subsec- 
tion IIIC. To capture the entire parameter range of the 
model, we investigate the properties of the effective spin 
Hamiltonian for positive and for negative exchange cou- 
pling. We start our investigations of the effective spin 
model in subsection HID by looking at its classical and 
semiclassical behavior. Surprisingly, we find that for pos- 
itive J there exists a very large manifold of degenerate 
classical ground states (GS). The semiclassical spin wave 
analysis (discussed in the subsection III E and limited to 
the most symmetric GSs) does favor some of those states, 
but does not give a definite answer concerning the real 
nature of the quantum ground state. In section IV we 
present the results of exact diagonalizations of finite cells 
of the realistic spin-1/2 version of our model. It turns 
out that even in this extreme quantum limit the ground 
state of our model exhibits long-range Neel order of the 
same structure as is found in the classical version. For 
positive exchange coupling, J > 0, we observe a very 
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high density of low-lying eigenstates of the effective spin 
model. We associate these low-lying states with the man- 
ifold of classical ground states whose degeneracy is lifted 
by quantum zero-point fluctuations. In the concluding 
section, Sec. V, we discuss experimental routes towards 
verification of our results, and detection of the predicted 
effects. The paper contains two appendices, in which we 
present the details of the calculations of the couplings 
in the Hubbard model, and the mean-field theory of the 
single component Bose gas, respectively. 

II. ATOMIC GASES IN KAGOME LATTICES 
A. Creation of optical kagome lattices 

In the following, we consider the atoms confined mag- 
netically or optically in the z direction at z = 0. The 
atoms form effectively a 2D system in an optical lattice 
in the x-y plane. In order to create a kagome lattice in 
this plane one can use red detuned lasers, so that the 
potential minima coincide with the laser intensity max- 
ima. A perfect triangular lattice can be easily created by 
two standing waves on the x-y plane, cos 2 (ki.2i"), with 
ki.2 = fc{l/2, ±\/3/2}, and an additional standing wave 
cos 2 (k 3 r + <j>), with k 3 = fc{0, 1}. The resulting trian- 
gles have a side of length 2ir/V3k. By varying (f> the 
third standing wave is shifted along the y axis, and, in 
principle, a kagome pattern could be realized. 

Unfortunately, this procedure presents two problems. 
First, three lasers on a plane cannot have mutually or- 
thogonal polarizations, and consequently undcsired in- 
terferences between different standing waves occur. This 
problem has, however, a relatively simple solution: un- 
desired interferences can be avoided by randomizing the 
relative orientation of the polarization between differ- 
ent standing waves, or by introducing small frequency 
mismatches, which, however, have to be larger than any 
other relevant frequencies. The second problem is much 
more serious, and is caused by the diffraction limit. Let 
us denote by £ the ratio between the separation between 
maxima of the laser intensity (i.e. minima of the resulting 
optical potential in the case of red detuned laser beam) 
and the half- width at half maximum (HWHM). To have 
a good resolution of the potential minima one needs £ to 
be definitely significantly larger than 2. In the case dis- 
cussed above, however, £ is only about 4 at <j> — ir in the 
ideal kagome case. Because of that, for any 0, the three 
potential minima forming the kagome triangles cannot 
be resolved. 

We propose to use the super-lattice technique [12] 
which we briefly describe in the following paragraphs, as 
a method to generate ideal and trimerized optical kagome 
lattices. The proposed experimental set-up is schemati- 
cally shown in Fig. I. There are three planes of standing- 
wave laser beams, and the wave vectors of these lasers lie 
on a plane. In the particular case of Fig. 1, we have three 
standing waves (a triple) in each plane. The laser fields 



within each plane are phase-locked. A kagome lattice will 
be formed by the intensity pattern that results from the 
sum of the laser intensities of the triples in the x-y plane. 

In order to resolve the three potential minima in the 
unit cell of the kagome lattice we must use at least two 
standing waves in each of the three vertical planes shown 
in Fig. I. While the wave-fields in the same plane must 
have identical polarizations, the fields in different planes 
should not interfere. As mentioned above, undesired in- 
terference cross-terms in the total intensity of the fields 
can be removed either by randomizing the relative orien- 
tations of the polarizations between waves from different 
planes, or by introducing small frequency mismatches. 
With this set-up consisting of 2 waves per vertical plane, 
we obtain the following intensity pattern in the x-y plane, 

3 

J(r) = I Q ^2 [cos(k,r + a^/2) 
i=i 

+ 2cos(k 4 r/3 + ( r 4 </>/6)] 2 , (r = (x,y)), (I) 

where ai = — 1 and o\ = cr 3 = I and the index i enu- 
merates the vertical planes. The pattern formed by the 
maxima of the intensity I(r) changes between a trian- 
gular lattice at <fi = 0, and trimerized kagome lattices 
with varying mesh width for < < 7r, until at = it 
the uniform kagome lattice is reached. In this limit one 
obtains the value £ « 7.6 at <j) — n. This is sufficient 
to create a well resolved ideal kagome lattice. Direct in- 
spection shows that in this case a moderately trimerized 
lattice can also be realized : £ remains sufficiently large 
for 57r/12 < (f> < 7r, so that the potential minima can still 
be resolved. 

With the additional third beam shown in Fig. I, a 
resulting intensity pattern 

3 

J(r) = h [ cos(k 4 r + 3<r^/2) 

+ 2 cos(k,r/3 + ct^/2) + 4 cos(k,r/9 + 2 , 

(r = (*,»)), (2) 

is obtained. With this arrangement it is possible to trans- 
form the optical potential smoothly from an ideal kagome 
case into a strongly trimerized lattice. The value of £ in- 
creases in this case to ~ 14, and remains large in wide 
range of angles (f>. 



B. Hubbard Hamiltonian 

Depending on the detuning of the laser relative to the 
resonance frequency of the atoms, either the minima or 
the maxima of the intensity patterns (1) and (2) form 
attractive potentials for the atoms. If these potentials 
are sufficiently strong, the tight binding approximation 
holds [47] , and the dynamics of the atomic gas can quite 
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Figure 1: (color online) Scheme of the proposed experimen- 
tal set-up. Each arrow depicts a wave vector of a standing 
wave laser. The three vertical planes intersect at an angle of 
120°. Dark (dark blue in the online version) spots in the right 
kagome figure indicate the potential lattice minima. 



generally be described by a Hubbard type Hamiltonian 
[2,3]: 

^Hubbard = - ^ kj{c\cj + c i c ]) + 7J ^ Unj(nj - 1) 

(y) i 
+ 2^2 U v n i n r ( 3 ) 

m 

Here c\ creates an atom in a Wannier state localized at 
the lattice site i. Depending on the atomic species the 
operators cj , Cj represent either fermionic or bosonic cre- 
ation and annihilation operators. The parameters iy, U 
and Uij of this Hamiltonian are matrix elements of the 
one-particle Hamiltonian and of the interaction poten- 
tials of the gas in the Wannier representation: 

Uj = (yVi\H \Wj), (4) 
where Hq is the one-particle Hamiltonian, 

F = -t^A + *,(!•), (5) 
2m 

with the one-particle potential v(r) cx ±J(r)- see (1) and 
(2). The sign depends on the detuning. For the Bose 
gas interacting via short range Van der Waals forces, the 
scattering at low energies occurs via the s-wave channel, 
and is adequately described by the zero-range potential, 
so that 

U = 5 2D / d 2 a;|W,(r)| 4 , (6) 

whereas 

Uij = 32D / d 2 x|W 4 (r)| 2 |W,(r)| 2 , (7) 

where the coupling (72D = / ^'Kh 2 a s /mW with m the 
atomic mass, and with W the effective transverse width 



of the 2D lattice in the z direction. In the case of polar- 
ized fermions U vanishes, since s-wavc scattering is not 
possible due to the Pauli principle. The nearest neigh- 
bors interaction, on the other hand, are possible, and in 
the case that they are due to dipolar forces (cf. [48]) or 
similar long range forces the couplings become 



Uij ~ y'd 2 a ;d 2 x'|W 4 (r)| 2 F(r-r')|W J (r')| 2 , (8) 



where V(r) is the interparticle potential. Obviously, the 
same expression holds also for bosons interacting via the 
potential V(r). The Hubbard Hamiltonian (3) does not 
necessarily describe the physics of bare particles; it may 
equally well describe the physics of composite objects, 
such as, for instance, composite fermions that arise in 
the analysis of Fermi-Bose mixtures in the lattice in the 
strong interaction limit [21]. The nearest neighbor inter- 
actions and tunnclings arc induced by the original hop- 
ping of bare fermions and bosons, and the corresponding 
values of tij and Uij have to be calculated from the bare 
couplings following the lines of Ref. [21]. 

In this paper we present explicit results for the tun- 
neling matrix elements and the interaction strengths 
U and Uij in the case of zero-range potential- expres- 
sions (4), (6), (7). To this aim we need to determine 
the Wannier functions Wj(r) for kagome type lattices. 
The method by which this task can be accomplished is 
presented in details in App. A. 

For the ideal kagome lattice we have successfully gen- 
erated the exact Wannier functions, and calculated the 
couplings accordingly. These results were then compared 
with the results of the variational method employing a 
Gaussian ansatz (for details see App. A). Fig. 2 compares 
the results calculated with Wannier functions and the 
Gaussian ansatz. For moderately strong potentials, say 
larger than two times the recoil energy E vec , the Gaus- 
sian approximation becomes appropriate, giving errors 
less than 50%. For sufficiently high potential amplitudes 
> 5E TCC , the results obtained with the Gaussian approxi- 
mation become practically indistinguishable from the ex- 
act Wannier results. 

Generating well localized Wannier functions in the 
trimerized lattices is a difficult task. For this reason, 
guided by the results for the ideal kagome lattice, we have 
limited ourselves here only to the results of the Gaussian 
approximation. Fig. 3 shows the hopping and interaction 
matrix elements depending on the trimerization angle. 
The perfect kagome lattice can be obtained by setting 
(j> = 7T. As expected, trimerization does not affect the 
on-site interactions very strongly, but does change the 
tunneling rates by orders of magnitude. Already a rela- 
tively moderate trimerization introduces large difference 
between the inter- and intra-trimer hopping elements. 
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Figure 2: (color online) Couplings for a perfect kagome lattice 
obtained using the Gaussian approximation (dotted lines) and 
using the exact Wannier functions (solid lines). Plot (a): hop- 
ping matrix elements, plot (b): contact interaction in units of 
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Figure 3: Couplings for a trimerized lattice obtained with a 
Gauss-function: (a) inter-trimer (</> < ir) and intra-trimer 
(4> > 7r) hopping matrix elements, (b) contact interaction 
terms in units of -E rC c- 



C. Bose gas in the trimerized kagome lattice 

As we have mentioned in Sec. I, in the present paper 
our main focus will be on a gas of spinless fermions on the 
trimerized kagome lattice at 2/3 filling. Other cases of 
interest which we will briefly discuss now include a Bose 
gas and a Fermi-Fermi mixture in this lattice. 

In order to facilitate the calculations, we add to the 
Hamiltonian (3) a term of the form — fJ,^2 i rii, where // is 
the chemical potential, that controls the average particle 
number of the system. Working with a fixed number of 
particles is possible, but technically very tedious. In the 
trimerized kagome lattice, the couplings tij take the val- 
ues i, t' for intra- and inter-trimer hopping, respectively. 
We set also f/y — V and = V , for intra- and inter-trimer 
interactions. 

In Ref. [11] we have considered the limiting case of 
hard core bosons, when U was much larger than any other 
energy scale, i.e. two bosons were not allowed at the same 
site. We have shown then that in the strongly trimerized 
case (t' , V <C V < t) the system will enter a trimerized 
Mott phase with the ground state corresponding to the 
product over (independent) trimers. Depending on the 
particular value of p = (/i — V)/(2t + V) we may have 
(p < -1), 1 (-1 < p < 0), 2 (0 < p < 1) or 3 (1 < p) 
bosons per trimer, i.e. filling factors v = 0, 1/3, 2/3 or 
1 boson per site. For fractional filling, the atoms within 
a trimer minimize the energy forming a, so-called, W- 



state [49]: \W) = (|001) + |010) + |100))/%/3 for v = 1/3, 
and \W) = QUO) + |101) + |011))/a/3 for v = 2/3. It is 
worth noticing that VF-states themselves have interesting 
applications for quantum information theory (c.f. [50]). 
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Figure 4: Mott phases (denoted by the corresponding particle 
numbers per trimer) of the state with lowest energy in the 
t — (fi + 2t) plane for zero inter-trimer hopping t' = 0. 




Figure 5: (color online) Phase boundaries between Mott- and 
superfluid-phase in parameter space of the hopping elements 
t, t' and the chemical potential fi. Below the loops the state 
is in a Mott-phase, where the number of bosons per trimer is 
displayed in the diagram. 

Generalizing the Landau mean- field theory of Ref. [19], 
we have obtained the phase diagram in the i' = t 1 /(2t + 
V) and p plane with characteristic lobes describing the 
boundaries of the Mott phases, given by F = (\p\ — l)/2 
for \fi\ > 1, and? = (3/2)|#|(1-|m|)/(4-|m|) for < 1. 
Observations of this Mott transition require tempera- 
tures T of the order of t' , i.e. smaller than t and V. 
Assuming that U is of the order of few recoil energies 
[3], this requires T to be in the range of tens of nK. The 
results for t < V are qualitatively similar. 

In this paper we present a method to generalize these 
results to the case, when the bosons are not necessar- 
ily hard core, i.e. U may be comparable with t. For 
simplicity we set EZy = 0, so that the Hamiltonian is 
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still described by the three parameters: t, t' for intra- 
and inter-trimer hopping, and U for the on-site inter- 
actions. Obviously, for vanishing inter-trimer hopping, 
t' = 0, the system is in a Mott-insulating state with a 
fixed number of particles per trimer. The corresponding 
Mott states are displayed in the phase diagram in the 
t—(fj, + 2t) plane in Fig. 4. As t' is increased the sys- 
tem undergoes a phase transition into a superfluid state. 
To obtain the phase diagram for this transition, Fig. 5, 
we have used a generalization of the Landau mean field 
approach of Fisher et al. [19, 51, 52], also investigated 
in [53]. Details of the method can be found in the Ap- 
pendix B. In our calculations we have confined ourselves 
to values of the chemical potential such that the particle 
number per trimer does not exceed four. In Fig. 5, further 
lobes with higher particle numbers will occur along the 
^i-axes for higher values of \x than those shown. Instead 
of calculating the so called super-fluid order parameter 
ip = (bi) = (b\) self-consistently a fully analytical expres- 
sions describing the boundaries in Fig. 5 can be obtained, 
as it is shown in the Appendix B. We mention that by 
using a cell strong coupling perturbative expansion [54] 
the phase boundary can be obtained with relatively little 
numerical effort with the accuracy of a Quantum Monte 
Carlo simulation. 



D. Fermi- Fermi mixture at 1/2 filling in the 
trimerized kagome lattice: Spin 1/2 Heisenberg 
ant ifer romagnet 

For a fermion-fermion mixture, instead of the extended 
Hubbard model described by the Hamiltonian (3) which 
includes nearest neighbor interactions U[j , it is more ap- 
propriate to consider the Hamiltonian with on-site inter- 
actions only: 

(ij) » 

Here fj and / denote the fcrmion annihilation opera- 
tors for the two species, and rij, hi are the correspond- 
ing occupation operators. The tunneling matrix elements 
are tij = t for intra-trimer and tij = t' for inter-trimer 
nearest-neighbor tunneling. Hpp is then the spin- 1/2 
Hubbard model. In the strong coupling limit, t, t' <C V, 
this model can be transformed into the t — J model [16] 
which reduces to the spin- 1/2 Heisenberg model for half 
filling, 

H FF -» H Heisenberg = J S i- S i 

{ij) intra 

(ij)inter 

where J = 4t 2 /V, J' = 4t ,2 /V, and S = (S x , S* S z ) 
with n-h = 2S Z , pf = S X + iS y , and pf = S x - iS y . 
It is exactly the model described by the }j Helsenber9 that 



has been studied by Mila and Mambrini [38, 39] in their 
effort to gain a physical understanding for the low-lying 
part of the spectrum of the kagome antifcrromagnet. The 
physics of this model is very interesting. In the trimerized 
case, it seems to be clear that the system qualifies as a 
RVB spin liquid of the second type. The large density of 
singlet and triplet excitations can be predicted quite well 
by analyzing the number of "relevant" dimer coverings 
of the trimerized lattice. The singlet-triplet gap, if it 
exists at all, is extremely small. All of these findings 
have so far no experimental confirmation. Experiments 
on this system are thus highly desirable, and we hope 
that ultracold atoms will allow to realize them. 



III. SPINLESS INTERACTING FERMI GAS IN 
THE TRIMERIZED KAGOME LATTICE AT 
2/3-FILLING 

A. Experimental realization 

Before we start to discuss the properties of the spin- 
less interacting Fermi gas in the trimerized kagome lat- 
tice, we shall first discuss the possibilities of preparing 
such a system. There are essentially two ways of achiev- 
ing this goal. First, we may consider an ultracold gas 
of fermions that interact via dipole-dipole forces. Bose- 
Einstcin condensation of a dipolar gas of chromium atoms 
has been recently achieved by the group of T. Pfau [55]. 
The (magnetic) dipolar interactions in chromium are sig- 
nificant, but not very strong. There are many ongoing 
experiments , however, aimed at the creation of ultracold 
gases of heteronuclear molecules, that could carry elec- 
tric dipole moments of the order of a Debye (cf. [56]). 
The observation of physics described in this paper using 
heteronuclear molecules with such strong dipolcs should 
be possible already at temperatures T ~ lOOnK. 

Another possibility of creating an interacting Fermi gas 
is to use the gas of composite fermions that appears in the 
low temperature behavior of Fermi-Bose mixtures in the 
limit of strong Bose-Bose and Bose-Fermi interactions. 
As we have mentioned above the physics of such compos- 
ite fermions is described also by an extended Hubbard 
model, in which the couplings result from virtual tun- 
neling processes involving bare fermions and bosons. In 
this case the observation of the low temperature physics 
requires achieving low, but not unrealistic temperatures 
T ~ 10 - 50nK (cf. Ref. [22]). 

The low energy states may be prepared by employ- 
ing adiabatic changes of the degree of trimerization of 
the lattice. For instance, one can start with a com- 
pletely trimerized lattice; then the filling v = 2/3 may 
be achieved by starting with v = 1 and by eliminating 
1 atom per trimer using, for instance, laser excitations. 
One can then increase t' and U 1 slowly on a time scale 
larger than 1/J (~seconds). Alternatively, one could 
start with v ~ 2/3 in the moderately trimerized regime. 
As in Ref. [3] , the inhomogeneity of the lattice due to the 
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trapping potential would then allow to achieve the Mott 
state with v = 2/3 per trimer in the center of the trap. 
Nearly perfect 2/3 filling can be reached by loading a 
BEC of molecules formed by 2 fermions into a triangular 
lattice, generating a MI state adiabatically, transforming 
the lattice to a trimerizcd kagome one, "dissociating" the 
molecules by changing the scattering length to negative 
values, and by finally optically pumping the atoms into 
a single internal state. Preparing v — 2/3 might involve 
undesired heating (due to optical pumping) which can 
be overcome by using laser or phonon cooling afterwards 
(cf. [57]). Note that the imperfections of v can be de- 
scribed by a "i — J" -kind of model, and are of interest 
themselves. 



B. Effective spin model 

The spinless Fermi gas in the trimerizcd kagome lattice 
is appropriately described by the Fermi- Hubbard Hamil- 
tonian 

Hfh = - ^(Ujfifj + h.c.) + ^2 UijUiiij - ^2 /Mi , 

(ij) (ij) » 

.(11) 

where tij and Uij take the values t and U for intra-trimer 
bonds and t' and U' inter-trimer bonds. \i is the chemical 
potential, and = are the occupation numbers with 
fi, fj the fermion annihilation and creation operators. In 
the following we denote the sites of each trimer by 1, 2, 
3 in the clockwise sense as shown in Fig. 6 




where n = p f and where n = n + + /!/- is the 

total number of fermions in the trimer. In the strongly 
trimerizcd limit the number of fermions is identical in 
each trimer. It is controlled by the chemical potential: 
for U + J < fi < 2U + J there are two particles in 
each trimer, one of them occupies the symmetric state, 
|1) = /*|0), while the second one occupies either one of 
the chiral states |1±) = 

In the inter-trimer part of the Hamiltonian 
Hhf, Eq. (11), we neglect the hopping term, 
~ T,( a i,/3j) t' (flj + h.c.) {a, /3 = 1, 2,3 refcring 
to intratrimer indices and i numbering trimers), since 
any real (first order) hopping process leads to an excited 
state whose energy is 0(U) higher than the ground- 
state energy and since second order (virtual) hopping 
processes yield contributions which are small of order 
t' 2 /U. Then, the inter-trimer part of Hhf reduces to 

jji 

ftinter _ — ^(ni^na^+Sj + ri,2,in 3 ,i+g 2 + n2,%ni,j + 6 3 
i 

n^in^i+Si + n 3ti n2,i+s 5 + ni,»n2,»+a 6 ) (I 3 ) 

Here, 8 U , v = 1, •• • ,6, denote the six vectors pointing 
from the central triangle to the six neighboring triangles, 
see Fig. 6. 

Next we express the occupation numbers n a ; , a — 
1, 2, 3, in terms of the fermion operators /, f± (we sup- 
press the site index): 



n 2 = 



n + (ft + /!)/ + fHU + /-) + r x ] , (14) 
n + (z+fl + z-fl)f + f\z-f+ + *+/_) 



+ cos(27t/3)t x + sin(27r/3)r ! ' 



(15) 



n 3 



1 



n + (4/t + zlfl_)f + P(z 2 _f + + 4/_) 



+ cos(27r/3)r x - sm(2n/3)T y ] . 



(16) 



Figure 6: The vectors Si, i = 1 . . . 6, pointing from the center 
of a given trimer to the centers of neighboring trimers. Num- 
bering of the sites of the trimers is shown in the triangle on 
the left. 

In this section it is our aim to derive from the Hamil- 
tonian (11) an effective spin Hamiltonian that captures 
in the strongly trimerizcd limit, t' ,U' <3jC t < U, the low- 
energy physics the model (11). 

The intra-trimer part of the Hamiltonian Hfh is diago- 
nalized by introducing instead of the local fermion modes 
fi, h and / 3 the symmetric mode / = (/!+ f 2 + /s)/\/3 
and the left and right chiral modes f± = (f\ + z±f 2 + 
4/3)/V3 with z± = cxp(i27r/3): 

Hfg a = -Mn + tn+^(n 2 - fi) - fin, (12) 



Here, the (pseudo-)spin operators t x :— |(/^/_ + 
/!/+), f» := /_ - /+/+) connect the right- and 
left-handed chiral fermion states. Inserting expressions 
(14) - (16) into H H lt F er , Eq. (13), yields bilinear terms 
in f x , t v , linear terms in f x and t v , bilinear terms in 
P, f and linear terms in p and /. Since none of these 
terms changes the total number of fermions in anyone of 
the trimers, we may set n = 2 in the resulting expression 
for H^p 7 '. However, terms containing the annihilation 
operator / promote the fermion in the symmetric state 
of a given trimer into the non-occupied chiral state of the 
same trimer. A glance at H l p H ra , Eq. (12), shows that 
the energy of this excited state is 0(t) above the ground- 
state energy. Thus, on account of analogous arguments 
as were given above for the neglect of the hopping term 
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of H l ftp , we also neglect all terms containing the oper- 
ators p , f. The linear terms in t?, rf sum to zero in 
the sum over the sites i so that we arrive at the following 
effective inter-trimer Hamiltonian (we omit an irrelevant 
constant): 



(a) 



7 f\ 



(b) , 



(c) 



H, 



eff 



T N 6 
»=1 i/=l 



Figure 7: (a) Classical ground state configuration for J < 
hs u )- (I?) (config. A), (b), (c) Classical ground state configurations 
for J > (configuration B and ferromagnetic configuration). 



Here, i are the sites of a triangular lattice of N sites 
on which the trimers are located, J = 4U'/9, and the 
vectors d„, v = 1, ■ ■ ■ 6, are the same as in Fig. 6. 
In Eq. (17), Ti(4>) = cos(0)ff + sin(0)rf and (j> ii g 1 = 
4>i,8e = 0, (j>i,5 2 = (j)^s 3 = 27r /3,^ 4>i,d4. = <t>i,6 B = 
-27r/ 3 >_ = 0i,« 2 = ^ 27r A 0i,i 3 = 0i,« 4 = 

and 4> iiS5 = 4> h s 6 = 27r/3. 



C. Effective spin model: relation to kagome 
ant ifer romagnet 

At this point it seems appropriate to briefly discuss the 
connection between the effective Hamiltonian H e ff de- 
rived here as model for the dynamics of fermionic atoms 
on a trimerized kagome lattice and the model Hamilto- 
nian that has been derived by Subrahmanyam [37] and 
has later been employed by Mila and Mambrini [38, 39] to 
explain the origin of the high density of low- lying singlets 
of the Heisenberg antiferromagnet (AF) on the kagome 
lattice. Mila considers the spin 1/2 Heisenberg model on 
the trimerized kagome lattice with a strong intra trimer 
coupling J and a weak inter trimer coupling J'. In the 
lowest order perturbation expansion with respect to J' 
he arrives at the effective Hamiltonian 



J' 



Hlrip- kag = ^ HijiSv)Hij{Th 



(18) 



(ij) 



where Hij(Sxj) = S V iS V3 and where Hij(r) is that mem- 
ber of our model H e f f that is associated with the bond ij. 
The operator S V i acts on the total spin of the trimer at 
site i, the trimers form a triangular lattice. In the deriva- 
tion Hlf } n - ka9 the Hilbert space of the three 5=1/2 
spins of the individual trimers has been restricted to the 
subspace of total spin 1/2 states. The four states of this 
subspace can be be specified by the z - component of their 
total spin and by two (spin)-chiralities. The Heisenberg 
type Hamiltonian Hij(S v ) acts on the two spin states of 
the trimers at sites i and j , Hij (r) acts on their chirali- 
ties. Obviously, ff^yp^ ka 9 turns into our model Hamil- 
tonian H e ff, if the trimer spins are fully polarized, 
e g Sf, = 1/2 for alH. This state can be reached by ap- 
plying a sufficiently strong magnetic field to the original 
trimerized kagome AF such that the total magnetization 
reaches 1/3 of the saturation magnetization, i. e. a mag- 
netic field that establishes the 1/3 magnetization plateau. 



quettes (19). 



D. Effective spin model: Classical aspects 

As is obvious from the derivation of the Hamiltonian 
(17), only its r = 1/2 quantum version can serve as a 
realistic effective model for the atomic Fermi gas in the 
trimerized kagome lattice. Nevertheless, for orientational 
purposes it is useful to first consider this model in the 
classical limit and to also calculate its excitation spec- 
trum in the semiclassical approximation, i. e. in the lin- 
ear spin-wave (LSW) approximation. We first describe 
the symmetries of the model Eq. (17). 

We have found that this model, is not only transla- 
tionally invariant, but is also invariant under the point 
group of order 6, Z 6 = Z 3 - Z 2 , where the generator of Z 3 
(order 3) is the combined rotation of the lattice by the 
angle 47r/3 and of the spins by the angle 27r/3 around 
the z axis, while the generator of Z 2 (order 2) is the spin 
inversion in the x-y plane, ff — ► —Tf, tJ l — > —ff. The 
model possesses no continuous spin rotational symmetry 
and the lines bisecting the angle between two adjacent 
lattice directions of the triangular lattice are not mirror 
lines. 

In Figs. 7a, b, c we show the three ordered classical 
states with small unit cells on the triangular lattice that 
are compatible with this point group symmetry: a ferro- 
magnetic state and two 120° Neel states labeled A and B 
which differ by the distributions of the chiralities x over 
the cells of the lattice as indicated by "+" and "-" signs. 
For an elementary cell of the triangular lattice whose cor- 
ners are labeled i,j,k in the counterclockwise sense, \ is 
defined as [58] 

Xijk = (rfT/-rfr;) + ( i,j^j,k ) + ( j,k^k,i ) 

(19) 

X is positive (negative) if the spin turns in the counter- 
clockwise (clockwise) sense as one moves around a trian- 
gular cell in the counter-clockwise sense. Because of the 
lack of mirror symmetry mentioned above it is not sur- 
prising that the two Neel states have different energies: 
^class = lr 2 JN, E^ s = - jT 2 JN. Here, N is the num- 
ber of sites and the superscripts A and B correspond to 
the labels of the Neel states in Figs. 7a, b. 

More surprisingly the ferromagnetic state is found to 
be degenerate with the Neel state B in the classical limit, 



^ciass = -^"dass- Furthermore, as is indicated by the an- 
gle 9 in Figs. 7a, b, the classical energies of the three 
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triangular lattice (see Fig. 7c): 



(C) 



/ / / / 



Figure 8: (color online) N=12 spin cell: (a) Configuration B; 
(b) Localized defect in configuration B marked by the trian- 
gular contour, (c) Ferromagnetic configuration; (d) Line de- 
fect in ferromagnetic configuration. The open arrows present 
spins that do not belong to the 12-spin cell. Their orientations 
are determined by the boundary conditions of the cell. 



structures do not depend on their direction relative to 
the lattice directions. In summary, in the classical ap- 
proximation the Neel state A is the ground state (GS) of 
model Eq. (17) for negative coupling, J < 0, while for 
positive J there are at least two classically degenerate 
GSs, the Neel state B and the ferromagnetic state. 

In fact, we have performed a numerical analysis of the 
12-spin cell by fixing the direction of every spin to nir/3 
with n = • • • 5, so that there were 6 12 classical spin con- 
figurations. This analysis has revealed that for J < 
there are 6 ground states each of them exhibiting the 
Neel order of type A (the six fold degeneracy comes from 
a Zq symmetry of our model). The results are dramat- 
ically different in the J > case, where we have found 
in total 240 degenerate classical GSs, among which the 
pure Neel states of type B and ferromagnetic states sum 
up to a small fraction. For an illustration, see Fig. 8 
where two ordered GS with very large unit cells (Figs. 
8b, 8d) together with their parent states (Figs. 8a, 8c) 
are presented. As will be seen below, the large number of 
degenerate classical GSs may find its analogue in a large 
density of low-lying excitations of the quantum version 
ofEq. (17). 



E. Effective spin model: spin wave theory 

The linear spin-wave expansion around the ferromag- 
netic GS based on the Holstein-Primakoff expansion [16] 
of the spin operators t x , f v is straightforward. The spin- 
wave frequency depends on the direction 9 of the mag- 
netization relative to the main lattice directions of the 



fcrro 



where 



,2tt 



/(q, 9) = — sin(— + 9) sin0 cos(q<5 1 ) 
o 

27T 2tt 
+ sin(y + 9) sin( — - 9) cos(q<5 2 ) 

2ir 

+ sin( 9) sin 9 cos(qJ 3 ) . 

3 

The quantum correction to the GS energy, SE fc "'°, is 
found as 

8E ic "°{9) = -5 J T + -L ^ w fero (q, 9). (20) 



(In the above expression and till the end of this sec- 
tion all expressions for energies present energies per site.) 
Evaluating the sum in this expression one finds that 
SE ie "° is minimal, if 9 takes one of the six values 7m/3, 
n = 0, • • • , 5. Thus, owing to the lowest order quantum 
corrections to the GS energy the magnetization of the 
ferromagnetic state locks in on one of the directions of 
the triangular lattice, i. e. the ferromagnetic GS be- 
comes sixfold degenerate in accordance with the order of 
the point group of our model, Eq. (17). With inclusion 
of these lowest order quantum corrections the GS energy 
for the preferred values of 9 is given by 



E tem = ~J[t(t- 
4 



1) 



0.901rl 



(21) 



In order to obtain the spin-wave frequencies for the two 
120° structures we closely follow the method devised by 
Jolicceur and Le Guillou [59] for the semiclassical treat- 
ment of the Heisenberg AF on the triangular lattice. 
Since the unit cells of both Neel structures A and B con- 
tain three sites one obtains in both cases three surfaces 
of spin wave frequencies in the magnetic Brillouin zone 
(BZ), and wf, a = 1, 2, 3 . As for the ferromagnetic 
state, the values of these frequencies depend on the an- 
gle 9 between these structures and the main directions 
of the triangular lattice, and hence the quantum correc- 
tions to the GS energies of the states A and B, SE A 
and SE A , depend on 9. For general 9 the expressions for 
the spin-wave frequencies are rather complicated. How- 
ever, by considering small deviations of 9 from the the 
values 7m/3 we find that as in the ferromagnetic case 
the quantum corrections SE A (9), 5E B {9) are minimal 
for 9 — irn/3, n = 0...5, i. e. the Neel structures A 
and B also lock in on the directions of the triangular lat- 
tice and hence both Neel states are sixfold degenerate. 
Remarkably, for the structure B all three branches of 
spin-wave frequencies are dispersionless with the lowest 
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branch consisting of N/3 zero modes, u>f (q) = for all q 
in the magnetic BZ. This is reminiscent of the Heisenberg 
model on the kagome lattice (HAK) for which one also 
finds N/3 zero-frequency spin wave modes. The nature 
of these zero modes is, however, quite different for the 
two models. While they correspond to simultaneous out- 
of-planc rotations of six-spin clusters in the HAK [60], 
they represent rigid in-plane rotations of the three spins 
on the corners of an elementary triangle in our model 
(17). 

Since it is of interest, we also note here the expression for 
the GS energy of the state B after the lowest quantum 
correction has been included: 

E B =-|j[t(t + 1)-1.48t] (22) 

Comparison of Eqs. (21) and (22) shows that quantum 
fluctuations lift the degeneracy of the purely classical 
states. In this scmiclassical approach it appears that 
the ferromagnetic state is the GS. We recall, however, 
that there is a very large manifold of classical GSs. In 
this manifold there may well be states which have lower 
energies than the two states that we have compared here. 



IV. NUMERICAL RESULTS 
A. Numerical method 

To describe the physics of spinless fermions on a 
trimcrized optical kagome lattice at filling 2/3 we need to 
consider the model (17) for spin r = 1/2, i. e. in the ex- 
treme quantum limit. Questions to be answered for this 
case are: (i)Is the GS of the model (17) an ordered state 
or is it a spin liquid either of type /, i. e. a state with- 
out broken symmetry, with exponentially fast decaying 
spin-pair correlations and a gap to the first excitation, or 
of type II, i. e. a kagome-like GS, again without bro- 
ken symmetry, with extremely short ranged correlations, 
but with a dense spectrum of excitations adjacent to the 
GS. (ii) What are the thermal properties of our system? 
After all, the model can only be realized at finite, albeit 
low temperatures. 

To find answers to these questions we have performed 
exact diagonalizations (ED) of the the Hamiltonian (17) 
for cells of N = 12, 15, 18, 21 and 24 sites using ARPACK 
routines [61]. The sizes of systems that can be studied 
by ED are restricted by the amount of memory space 
that is required for storing the non-zero matrix elements. 
To reduce this requirement we block-diagonalized the 
Hamiltonian (17) by exploiting its invariance under TV- 
fold translations. It allowed us to reduce the problem of 
diagonalization of 2 N x 2 N matrix to N independent di- 
agonalizations of matrices of size ~ 2 N /N x 2 N /N. This 
simplification not only lowers the memory requirements 
but also greatly reduces the time of calculation, especially 
when a large number of excited eigenstates is of interest. 



Nevertheless, ED of this Hamiltonian remains a de- 
manding task, as in contrast to SU(2) invariant spin 
models the Hilbert space of the Hamiltonian (17) can- 
not be separated into subspaces of states with fixed total 
spin and total z-component of the spin. Because of this 
last circumstance we had to limit our study to systems 
of at most 24 spins. Fortunately, our results for 21 and 
24 spins show qualitative and quantitative resemblance. 
Therefore we regard them as representative for larger 
systems too. In presenting our results we shall mainly 
confine ourselves to the two largest systems, since the 
results for smaller systems suffer from strong finite size 
effects. We remark that only the 12- and the 21-site cell 
can be chosen such that these systems possess the full 
point group symmetry of the infinite lattice. The lack of 
this symmetry for the 15- and the 18-site cell adds to the 
large finite size effects observed in the results for these 
cells. 



B. Ground state and low temperature properties 

For J < 0, i. e. for attractive interaction U' between 
fermions on nearest neighbor trimers, the highest-levels 
of H e f f I J and the corresponding eigenstates are physi- 
cally most relevant. As will be seen below theses levels 
are well separated from each other so that we only need 
to calculate a few of them. The situation is drastically 
different in the case J > 0, where we need the low-lying 
states of H e ff. It turns out that there is an abundance 
of such low lying states. In this respect the spectrum of 
H e ff is reminiscent of the spectrum of the Heisenberg 
Hamiltonian on the kagome lattice [40, 42]. The answer 
to the question of whether there is long range order in 
our model (17) is found in Tables I and II, where we 
show our numerical results for the spatial spin-spin cor- 
relations for the J < and J > cases, respectively. 
The cells to which these tables refer are shown in Figs. 
9a, b, c. 

Let us first consider the case J < 0, Table I. We have 
not done a systematic finite-size analysis for these corre- 
lations. However, comparing the data for the quantum 
t = 1/2 systems with the classical correlations there can 
be little doubt that in its GS the system orders in the 
planar 120° Neel structure. The smallness of the out-of 
plane correlations lends further support to this conclu- 
sion. We have also calculated the expectation values of 
the chirality Xijk, Eq. (19), in the GSs of the 12- and of 
the 21-site cell and have found perfect agreement with the 
pattern of positive and negative chiralities of the classical 
configuration, Fig. 7a. Apparently, for J < quantum 
fluctuations have a rather weak effect on the GS proper- 
ties of our model (17). 

For the case J > our results for the spin-spin correla- 
tions are presented in Table II. As for the case J < we 
have not been able to perform a finite size analysis but 
again we interpret the data in Table II as evidence for 
the existence of planar 120° Neel order in the GS of our 
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(a) (b) 




Figure 9: (a) 12 spin cell, (b) 21 spin cell, (c) 24spin cell. Oj 
V, and □ mark the three sublattices. Primed sites belong to 
periodic repetitions of the cell containing the unprimed sites. 





1st neighbs. 
d = 1 


2nd neighbs. 

d = V3 


3rd neighbs. 
d = 2 


4th neighbs. 
d = V7 


classical 


-0.125 


0.25 


-0.125 


-0.125 


N = 12 


-0.137 


0.251 


-0.125 




TV = 21 


-0.134 
(-0.029) 


0.237 
(-0.004) 


-0.117 
(-0.004) 


-0.116 
(-0.003) 



Table I: Spin-spin correlations, (tq + Tq tJ) for J < . In 
the last row the r z — r z correlations are also shown (numbers in 
parentheses). Owing to the Zq symmetry of the Hamiltonian 
(17) the correlations depend only on the distance d from the 
central site (see Fig. 9). 



model (17). This contradicts the prediction of the LSW 
analysis according to which one might have expected to 
find a ferromagnetic GS (see subsection HIE). However, 
one must recall that besides the 120° Neel GS and the fer- 
romagnetic GS there are many more classical GSs. In a 
complete LSW analysis one would have to consider every 
one of these states, a task that is practically impossi- 
ble to perform. Relative to the in-plane-correlations the 
magnitude of out-of-plane correlations, which are not dis- 
played in Table II, is even smaller here than in the case 
J < 0. Further support for long-range order in the GS of 
the model (17) comes from a comparison of the spin-spin 
correlations of this model with the same correlations of 
the = 1/2 Heisenberg AF on the triangular lattice (TAF) 
which we have included in Table II. It is seen there that 
the GS correlations of the model (17) decay more slowly 
than those of the GS of the TAF which is known to pos- 
sess long range 120° Neel order [33]. 

Additional strong support for existence of a Neel or- 
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Figure 10: (color online) The curves from a bottom to top 
correspond to in-plane spin-spin correlation to the 1st neigh- 
bor (black, dash-dotted), 3rd neighbor (green, dashed), 4th 
neighbor (blue, dotted) and to the 2nd neighbor (red, solid). 
Two thousand lowest eigenstates were used in this calculation. 
The system size is JV = 21. 



dered GS in both J > and J < cases comes from an 
investigation of chirality patterns. In both situations the 
quantum mechanical calculation reveals that there exists 
a perfectly periodic pattern of chiralities Xijk as m the 
classical result. For J > we found that Xijk ~ TO-5 
while for J < Xijk ~ ±0.69. Both results are obtained 
in N = 21, where the Xijk was calculated for six trian- 
gles located around the central site (Fig. 7). The =F, 
± notation indicates opposite chiralities between J > 
and J < results. A comparison of these values to =f0.65 
(Neel B configuration) and ±0.65 (Neel A configuration), 
leaves little doubt on the nature of these GSs. Finally, 
please notice excellent agreement between quantum and 
classical calculation for J < 0. 

Values for the spin-spin correlations of the model (17) 
for finite albeit small temperatures are also displayed in 
Table II. For T = 0.005J about 800 low lying eigen- 
states were needed to achieve convergence in the data 
for the correlations. Although these finite temperature 
correlations are smaller in magnitude than the GS cor- 
relations, they decay as slowly with the distance as the 
GS correlations i. e. long-range order persists at finite 
temperatures. This is not surprising since our model (17) 
has no continuous symmetry. One thus expects the order 
to vanish at a finite temperature T c in a first or second 
order phase transition. The temperature dependence of 
spin-spin correlations in A = 21 system is depicted in 
Fig. 10. 

The finite size effects affect the correlations very 
strongly for system sizes N < 21. In Fig. 12 we plot 
the spin-spin correlations for the various system sizes. 
The data for N=15, 18 are averages of the raw data for 
fixed lattice distances over the lattice directions. Be- 
cause of boundary effects the correlations do not show 
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Figure 11: Ground state energies in units of |J| as a function 
of the system size N. Solid (dashed) line is a linear fit to J > 
(J < 0) data. The fit gives in J > case: —0.22 J N — 0.07 J, 
while in the J < case: 0.40 J N + 0.20 J. 
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Figure 12: (color online) Spin-spin correlations for the various 
system sizes. The data for N=15, 18 are averages of the raw 
data for fixed lattice distances over the lattice directions. 



the expected six fold symmetry. As a consequence the 
data for N=15, 18 cannot be used in a finite-size ex- 
trapolation. Nevertheless, despite these large finite size 
effects, for both cases J > and J < the GS ener- 
gies can be reliably extracted from the data for all the 
cell sizes including the smaller ones. From the linear 
fits shown in Figs. 11a, b we obtain E& s = -0.40| J| as 
the GS energy in case A. This is to be compared with 
the classical GS energy (see subsection HID) 
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class 



-0.375|J|, (t = 1/2). In the same way we 
-0.22 J as the GS energy per site in case B, 
which is to be compared to the classical GS energy (see 
subsection HID) Eg^ = -|r 2 J = -0.1875 J. 
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Figure 13: (color online) Accumulated density of states, 
i. e. the number of states in the energy interval AE above 
the ground state for (a) J < and (b) J > 0. The curves 
in plot (b): blue (solid) N=24, red (dotted) N=21 and green 
(dashed) N=18. 



C. Low energy spectra 

Let us finally discuss the energy spectra of our model 
for both cases J < (A), and J > (B). Figs. 13a, 
b show the accumulated density of states (accumulated 
DOS) of our model (17) for the two cases. On account of 
the breaking of the discrete symmetries of the Hamilto- 
nian (17) by the 120° Neel order the standard expectation 
would be that the GS is sixfold degenerate for the infi- 
nite model and, since there is no continuous symmetry 
that could be broken, the excitations should be separated 
from the GS by a finite gap of the order of J. For finite 
systems the GS degeneracy will be lifted. Nevertheless, 
we expected to find six low-lying states in the gap below 
the lowest excited state. Fig. 13a, J < 0, does not re- 
flect this scenario convincingly. However, there are only 
a few states with energies substantially below 0.5 1 J\. We 
take this as an indication of a gap of this order of mag- 
nitude in the spectrum of the Hamiltonian (17) in the 
thermodynamic limit N — > oo. 

Obviously, for J > the spectrum differs drastically 
from the above expectations, see Fig. 13b. There is an 
abundance of very low-lying excitations, e. g. for N = 21 
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Figure 14: Specific heat for the N = 21- Eq. (23). Two 
thousand lowest eigenstates were used in the formula (23) to 
prepare this plot. 

there are about 2000 (800) states with energies less than 
0.09J (0.05 J) above the ground state. From the perfect 
symmetry of a finite temperature spin-spin correlations 
and their relatively slow decay with temperature, we con- 
clude that majority of these excited eigenstates support 
the spin order of the GS. 

Comparison of the lower panel of Fig. 13b with the 
scenario outlined above suggests that the gap, if any, is 
smaller than 0.5 • 10~ 2 J. 

The rapid increase of the accumulated DOS that sets 
in at excitation energies of this order of magnitude leads 
to peaks in the specific heat, 

1 d E t ^cxp[-^/(fcT)] 

NOT £.exp[-iV(fcT)] ' 1 1 

at the corresponding temperature. We have checked that 
the peak shifts towards lower temperatures as the size of 
the system increases. Indeed, in the N — 24 system, we 
found the peak at kT w 2.5 • 10~ 3 J while for N = 21 
it is at kT w 3.6 • 10~ 3 J, see Fig. 14 for the N = 21 
results. The precise determination of the peak position 
and amplitude in the N = 24 system requires, however, 
a calculation based on more excited eigenstates than we 
have been able to get (~ 240, see Fig. 13b). The shift of 
the peak position between these two systems reflects the 
slightly different behavior of their accumulated DOS. It 
would be very interesting to know, whether this trend 
continues for larger N so that in the thermodynamic 
limit, N — > oo, the specific heat no longer decreases to 
zero at T = 0. In this case our model (17) would be an ex- 
ample of a quantum model with a finite zero-temperature 
entropy. 

Because of the strong finite size effects in the data for 
N < 21 the growth law of the accumulated DOS with 
the system size N cannot be extracted reliably from our 
data. They are compatible, however, with an exponential 
increase of the number of low-lying states with N. 

The features of the low-energy part of the spectrum of 
our model, Eq. (17), are strongly reminiscent of the low- 
energy part of the spin- 1/2 Hcisenberg antiferromagnet 



on the kagome lattice (HAK) [40, 42]. There is, however, 
one decisive difference between the two models: while all 
GS correlations were found to be extremely short ranged 
[43] in the HAK there is, in all probability, long range 
spin order in the GS of our model. The absence of long 
range order in the HAK led Mila and Mambrini [38, 39] 
to study the trimerized HAK in the basis consisting of 
all independent dimer coverings of the lattice by exclu- 
sively nearest-neighbor singlet pairs. By definition this 
restricted basis cannot produce any long range order in 
the GS of the HAK. Using it in analytical and in numer- 
ical calculations Mila and Mambrini were able to repro- 
duce the low- lying part of the spectrum of the HAK. In 
particular, they were able to determine the constant a 
in the growth law a N that describes the increase of the 
number of low-lying states in the HAK. However, the ap- 
proach of Mila and Mambrini is not suited for the treat- 
ment of our model for at least two reasons: (i) in contrast 
to the HAK our model is not SU(2) invariant. Therefore 
a restriction of the full Hilbcrt space of the model to ex- 
clusively singlet states is unwarranted, (ii) We need to 
describe spin-ordered states, and this is not possible in 
a basis consisting of products of nearest-neighbor singlet 
pairs. We suggest that for our model (17) the abun- 
dance of low-lying quantum states corresponds to the 
abundance of classical GS described in subsection HID. 
Zero-point fluctuations lift the degeneracy of the clas- 
sical states leaving the spin correlations that are built 
into these classical states qualitatively untouched. On 
account of its low-energy properties we have proposed 
the name quantum spin-liquid crystal for our system. 



V. CONCLUSIONS 

In this paper we have discussed in detail the physics of 
ultracold gases in trimerized kagome lattices. Observa- 
tion of this kind of physics, and detection of the predicted 
effects requires various steps: preparation of the trimer- 
ized lattice, loading of the considered gases, and detec- 
tion. The first step, i.e. the preparation of the kagome 
lattice, is discussed in detail in subsection II A. 

Probably, the easiest experiment to perform concerns 
the observation of the novel Mott phases of the Bose gas. 
Temperature requirements (~ lOOnK) are rather mod- 
erate. The challenging problem here is how to achieve 
1/3, 2/3 fillings, etc. In principle physics should do it 
for us, since the "exotic" Mott phases are the thermody- 
namic phases of the system at zero temperature. There 
is, however, another elegant method of preparing such 
phases. To this aim one should start with a triangular 
lattice and achieve a Mott state with 1, 2, 3, ... atoms 
per site. Then one should deform the lattice to a trimer- 
ized kagome. The detection of such Mott phases can be 
done simply by releasing the atoms from the lattice, as 
in Ref. [3]. Coherence on the trimer level will then be 
visible in the appearance of interference fringes in time- 
of-flight images, which should reflect the on-trimer mo- 
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Table II: Spin-spin correlations as in Table I, but for J > 
0. Sites j are numbered as in Fig. 9. For comparison the 
correlations (SESf + S%S]} for the spin- 1/2 Heisenberg AF 
on the triangular lattice are also displayed (data from [41]). 
Dots below a value for the correlation indicate that this value 
occurs repeatedly. 



mentum distribution <~ . cos k(fi — fj), where fj are 
the positions of the minima in a trimer. In spite of the 
appearance of these fringes, the Mott-insulator nature of 
the state would be apparent in the presence of a gap for 
the excitations, which can be observed by tilting experi- 
ments as those of Ref. [3] . The opening of the gap should 
be analyzed as a function of the trimerization degree t'/t, 
which can be controlled as discussed in Sec. II A. 

The experiment with Fermi-Fermi mixture is more de- 
manding. The main problem is, of course, the prepara- 
tion of the states in the low energy singlet sector. One 
possible way to prepare a singlet state in the trimerized 
kagome lattice with T < 3i/4 could employ the recently 
obtained Bose-Einstein condensates of molecules consist- 
ing of two fermionic atoms [62] at temperatures of the 



order of 10 nK. Such BECs should be loaded onto an 
ideal and weak kagome lattice. Note that the molecules 
formed after sweeping across a Feshbach resonance, are 
in a singlet state of the pseudo-spin s. This can easily be 
seen, because the two fermions enter the resonance from 
the s-wave scattering channel (i.e. in the symmetric state 
with respect to the spatial coordinates), and thus are in a 
singlet state of the pseudo-spin (i.e. antisymmetric state 
with respect to exchange of electronic and nuclear spins) . 
Since the interaction leading to the spin flipping at the 
Feshbach resonance [63] is symmetric under the simul- 
taneous interchange of both nuclear and electronic spin, 
then the formed molecule remains in a pseudo-spin sin- 
glet. The typical size of the molecule is of the order of the 
s-wave scattering length a, and thus can be modified at 
the resonance [64], being chosen comparable to the lat- 
tice wavelength. Growing the lattice breaks the molecule 
into two separate fermionic atoms in neighboring sites in 
the singlet pseudo-spin state. In this way, a singlet cov- 
ering of the kagome lattice may be achieved, allowing for 
the direct generation of a RVB state [65]. 

Assuming that we can prepare the system in a sin- 
glet state at J' < T < J, then the density of states of 
the low lying singlet levels can be obtained by repeated 
measurements of the system energy. The latter can be 
achieved by simply releasing the lattice, so that after tak- 
ing care of the zero point energy, all of the interaction en- 
ergy is transformed into kinetic energy. In a similar way 
we can measure the mean value and the distribution of 
any nearest-neighbor two-spin correlation functions. To 
this aim one has to apply at the moment of the trap re- 
lease a chosen nearest neighbor two-spin Hamiltonian and 
keep it acting during the cloud expansion (for details see 
[25]). In a similar manner we can measure the spectrum 
of triplet excitation, by exciting a triplet state, which 
can be done by flipping one spin using a combination of 
superlattice methods and laser excitation [66] . The mea- 
surement of the singlet-triplet gap requires a resolution 
better than J'. 

A similar type of measurements can be performed in 
the ideal kagome lattice, when J = J 1 . In this case, 
the singlet-triplet gap is filled with singlet excitations 
[42] . By varying <f>, one can transform adiabatically from 
strongly trimerized to ideal Kagome, for which the final 
value of J will be smaller than the initial J, but larger 
than the initial J' . In that case, the system should re- 
main within the lowest set of 1.15^ states that originally 
formed the lowest singlet band. The singlet-triplet gap, if 
any, is estimated to be < J/20, and should be measurable 
using the methods described above. 

The observation of properties of the spinless interact- 
ing Fermi gas is also experimentally very challenging. 
The first step is to create the interacting Fermi gas, 
obviously. As we discussed it in subsection III A this 
can be achieved either with dipolar particles, or compos- 
ite fermions. Both of these possibilities are challenging 
themselves, although the rapid progress in cooling and 
trapping of dipolar atoms [55] and molecules allows one 
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to hope that interacting spinless Fermi gases will be rou- 
tinely available in the next future. Preparing of the 2/3 
filling is also a challenge, but several routes have been 
proposed in subsection III A. Yet another challenge is to 
measure the predicted properties of the quantum spin- 
liquid crystal. 

One quantity which should be possible to measure rel- 
atively easy, is the energy of the system. This can be 
done simply by opening the lattice; by repeated mea- 
surement of the energy E(T) at (definite) finite temper- 
atures one would get in this way an access to the den- 
sity of modes, i.e. one could compare the results with 
Fig. 13. From such measurements one could infer the 
existence of a gap -Egap, since, if E gap is large enough, 
E(T) becomes T- independent for kT < E gap . Various 
other correlations could be measured using the methods 
proposed in Ref. [25]. In order to measure planar spin 
correlations, one has, however, to lift the degeneracy of 
the f± modes, e. g. by slightly modifying the intensity 
of one of the superlattices forming the trimerized lattice. 
This should be done on a time scale faster than the char- 
acteristic time scales of other interactions, so that the 
state of the system would not change during the mea- 
surement. In such a case one can use far off resonant 
Raman scattering (or scattering of matter waves) to mea- 
sure the dynamic structure factor, which is proportional 
to the spatio-temporal Fourier transform of the density- 
density correlations. At frequencies close to the two pho- 
ton Raman resonance between the f± modes, only /+-/_ 
transitions contribute to the signal, and hence such mea- 
surement yields the desired information about the cor- 
relations ^ f+ an d the spin correlations of 
Figs. 10,12. 
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reads H = —t{c\c2 + c3,C3 + cXxi\ + h.c}. The eigenfunc- 
tions of H arc a = (1, 1, 1) T /V3 , oi = (-2, 1, 1) T /2 
and a 2 = (0, — 1, 1) T / y/2 with respective eigenvalues 
{—2t,t,t}. The transformation matrix (So, 01,02) 1 
leads to the states for which just one site of the trian- 
gle is occupied. A similar scheme can also be applied for 
the Bloch functions V'l( r ) = c ~ lkr u k (r), where k is the 
quasi- momentum, and u k (r) are the periodic functions 
of band i. For a particular quasi momentum k, the val- 
ues of w k (r) at the potential minima are ay(k), where 
j € {1,2,3} denotes the minimum within the trimcr 
(Fig. 6). Inverting the matrix for a particular k one 
obtains complex coefficients c op , which are then used to 
construct a periodic function having its maximum in only 
one of the three potential minima: w£ = E, c iU k , where 
o = {1, 2, 3} denotes the corners of the trimer. Similarly 
as the u k , the w k are functions with the same period- 
icity as the lattice. Summation of these functions over 
k with a proper phase leads to the Wannier functions 
[67]: = ^^k eikR,i, k. whcrc R denotes the posi- 

tion of the particular trimcr on which the maximum of 
the Wannier functions is located. 

The point group symmetry of the lattice is broken due 
to the choice of the particular set of basis vectors for the 
reciprocal lattice. A direct consequence of this fact is that 
the Wannier functions within a trimer cannot be trans- 
formed into each other by a rotation of ±2/37r around 
the center of the trimer. Hence one obtains different hop- 
ping probabilities between the sites of the triangle, tij — 
(Wk\H \Wi) = ^ Ek E M where e£ is the en- 

ergy for the quasi momentum k in band ji. The Wan- 
nier functions can be symmetrized by summing up Bloch 
functions multiplied with a k-dependent phase factor 
c~ lkri , where is the position of one of the three poten- 
tial minima within a cell. The hopping elements change 
then to = (I/AT, 2 ) E k cos k(ri - E G c i, Cj, e k , 
which are independent of the position now. The cost of 
the symmetrization is that the Wannier functions are not 
orthogonal anymore, but the overlap remains relatively 
small. 



Appendix A: CALCULATION OF THE 
PARAMETERS FOR THE HUBBARD 
HAMILTONIAN 

1. Wannier functions 

For periodic boundary conditions the linear part of the 
Hamiltonian (5) can be diagonalized in the quasi momen- 
tum space using the scheme of Bloch [47] . In the kagome 
lattice a single cell contains three equivalent potential 
minima, and hence three Wannier functions per unitary 
cell are required, which can be obtained by transforming 
the Bloch states of the first three bands into the Wannier 
basis. Let us first consider one particle that is placed 
in an isolated trimer. The Hamiltonian for this model 



2. Gaussian ansatz 

Apart from the case of kagome lattice, it is difficult 
to obtain the Wannier functions reliably. The coeffi- 
cients for the Hubbard-Hamiltonian can be alternatively 
obtained using a Gaussian Ansatz [68], which, in the 
case of a perfect kagome lattice, and for deep lattice po- 
tentials > 5E rcc leads to results which are practically 
indistinguishable from those of the Wannier functions. 
The shape of the Gaussian function reads f(x,y) — 
y/2/(cr x <T y n) exp (-x 2 /cr 2 ) exp {-(y - yo) 2 /^). The 
center yo, the widths o x and cr y are the varia- 
tional parameters minimizing the energy functional: 

E = /r co dx/r oo dy[V(/(x,y)) 2 + /( a ;,y) 2 V-(^y)]- 
The Gauss functions for the other two minima in the 
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trimer are obtained by rotating the Gaussian function 
by ±=7r around the center of the trimer. 



Appendix B: MEAN-FIELD THEORY FOR A 
BOSONIC GAS 



The boundaries between Mott-insulator and superfiuid 
phases can be obtained by means of a mean-field ap- 
proach similar as that employed in Ref. [21]. We con- 
sider only on-site contact-interaction terms, but contrary 
to Ref. [21] we do not restrict ourselves to the hard-core 
limit. The system is governed by a Bosc- Hubbard Hamil- 
tonian of the form: H = H tr + Hh op , with 



H tr = = -t^y {b\b, + h.c.) 



H hop = -t'Y(blbi + h.c), 



(kl) 



(Bl) 



(B2) 



where t, t' (denoting the intra- and inter-trimer hop- 
pings) and n are measured in units of the on-site inter- 
action potential U. 

Assuming a fixed number of atoms n per trimer, we 
consider all possible Fock-states of the form |nin 2 n 3 ) 
with m + n 2 + n 3 = n. For example, for one par- 
ticle per trimer the Hilbert-space contains the Fock- 
states 1 100) , |010) and |001). Since the model is in- 
variant under rotation of 2n/3 the eigenstates are of 
the form |Wi) = (|100) + z|010) + z 2 |001))/\/3, with 
z € { 1, exp (i§7r), exp (— i§7r)}, implying states with no, 
left- and right-chirality, also known as W states [49] . We 
denote z± = cxp(±i27r/3) and introduce the operators 

bs)/V3. 

a,(3 = 



B± = (6i + z ± b 2 + 4&s) /V3, B = {bi + b 2 



B a ,B^ 



= 5, 



a/3 



Their commutation relations are 
{0, +, — }. The chirality operator is defined as x = 
(B f + B+ - B^B-) mod 3. Eq. (Bl) can be rewritten 
into the form 



For a given number of particles n per trimer, we de- 
note by \W£) a particular normalized non-chiral sum of 
permutations of a set ni,ri2,n3, where fi denotes differ- 
ent allowed non chiral states. E.g. for two particles per 
trimer, two possible states are allowed: \W 2 ) = (|H0) + 
|101) + |011))/V3 and |TU 2 2 ) = (|200) + |020) + |002))/\/5. 
We can then diagonalize the Hamiltonian H tr in this ba- 
sis, H'jf = (W™\H tI \WP ), obtaining the eigenenergies e l n 
and eigenstates \i(> l n ), where < I < n. The lowest ener- 
gies e° for each particle number n have to be compared 
to obtain the ground-state in the (i-^)-phase space. 

If the inter-trimer hopping t' is small, the phase bound- 
aries in the t-t'-fi phase diagram can be well estimated 
by using a mean-field approach [19, 51, 52] We introduce 
the superfiuid order parameter ip = (bi) = (b\), for every 
site i. Neglecting fluctuations of bi, £>t in the second or- 
der, we can substitute bjbi — ip(bj + bi) — ^ 2 , and hence 
Hh op can be decomposed into a sum of single-site Hamil- 
tonians of the form: H hop w 6t'ip 2 - 2V3t'ip(B + B r ). 
The Hamiltonian H can be decomposed then into two 
parts H = H a + V, with H a = H tT + 6t'ip 2 and V = 
-2t'V3tp(B + Bl), where H is perturbed by V. Up to 
second order perturbation theory, the energy becomes of 
the form E = e° + rip 2 , where 



w^ 2 



E 

m— nil, 2 



(B4) 



The Mott-insulator to superfiuid transition may be iden- 
tified by the equation r = 0, since for r > the energy 
is minimized for ip 2 is zero, and for r < tp acquires a 
finite value. The equation r = defines a 2D manifold 
in the i'-i-/i-parameter space. 

As an example, we determine in this Appendix the 
expression for the boundaries of the Mott-phase with one 
particle per trimer. Due to the form of Eq. (B4) this 
calculation demands the knowledge of the eigenenergies 
and cigenfunctions for n = 1 and n = 2. For n = 1, 
= (|001) + |010) + |100>) /V3 = B+I^o), and e? = 
($>|fltr|V?) = -p-2t. Forn = 2, 



= i (l T V(l + 2t) 2 + 32t 2 ) -t-2/x, 
H tr = -tWlB + (t - it) [b\Bq + B\.B+ + and |^ 4 ) = cos0 o ,i|W 2 2 ) + sin fo^Wj), with 



(B5) 



1 -[(b£ + 2B\Bl) (B 2 +2B + B_) 
+ (s| 2 + 25^) (B\ + 2B Q B_) 
+ (Bl 2 + 2B ] q B\^ (B 2 _ + 2B B+) } . 



(B3) 



Therefore expi2n/3(B f + B + - B ] _B-),H tr ^ = 0, and 

hence the chirality of a state is a conserved quantity. It 
can be shown that the ground state has chirality zero, 
and therefore we restrict ourselves henceforth to these 
states. 



tan 0o.i = — ^{(1 + 2t) T y/(l + 2t) 2 + 32t 2 }. (B6) 

At t' = the region of 1 particle per trimer is provided 
by the condition < e° < e®, i.e. when — 2t < fj, < 
t + (1 - v/(l + 2t) 2 + 32i 2 )/2. 

After a straightforward but tedious calculation, we can 
then calculate the sum in Eq. (B4) 

m\v\rj\ 2 



E 

m=0,2,i 



(6/x - 24t - 4) 3 s 

- \xi2t + 1) - 8t 2 ~ 2t + yL / 



(B7) 
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Hence, solving for r = 0, we obtain the value of t' at the 
phase boundary: 

= 1/2 (^ 2 - M (2t+l)-8t 2 ) (2t + ») 

+ 8i)(2t + 1/3) - /z 2 - 8t 2 ' 1 ' 
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